Black hole entropy from entanglement: A review 
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Abstract 

We review aspects of the thermodynamics of black holes and in particular take into account the fact 
1—5 \ that the quantum entanglement between the degrees of freedom of a scalar field, traced inside the event 

■ horizon, can be the origin of black hole entropy. The main reason behind such a plausibility is that the 
well-known Bekenstein-Hawking entropy- area proportionality — the so-called 'area law' of black hole 

' q' physics — holds for entanglement entropy as well, provided the scalar field is in its ground state, or in 

other minimum uncertainty states, such as a generic coherent state or squeezed state. However, when 
1 ■ the field is either in an excited state or in a state which is a superposition of ground and excited states, 

(3JT)! a power-law correction to the area law is shown to exist. Such a correction term falls off with increasing 

area, so that eventually the area law is recovered for large enough horizon area. On ascertaining the 
I | location of the microscopic degrees of freedom that lead to the entanglement entropy of black holes, it 

■ is found that although the degrees of freedom close to the horizon contribute most to the total entropy, 
fS| \ the contributions from those that are far from the horizon are more significant for excited/superposed 

states than for the ground state. Thus, the deviations from the area law for excited/superposed states 
may, in a way, be attributed to the far-away degrees of freedom. Finally, taking the scalar field (which 
is traced over) to be massive, we explore the changes on the area law due to the mass. Although most 
of our computations are done in flat space-time with a hypothetical spherical region, considered to be 
the analogue of the horizon, we show that our results hold as well in curved space-times representing 
| static asymptotically flat spherical black holes with single horizon. 

' PACS Nos.: 04.60.-m, 04.70.-s, 04.70.Dy, 03.65.Ud 

x : 

1 Introduction 

One of the most remarkable features of black hole physics is the realization that black holes behave as 
thermodynamic systems and possess entropy and temperature. The pioneering works in the field of black 
hole thermodynamics started with Bekenstein [1] , who argued that the universal applicability of the second 
law of thermodynamics rests on the fact that a black hole must possess an entropy (S BH ) proportional 
to the area (An) of its horizon. The macroscopic properties of black holes were subsequently formalized 
by Bardeen, Carter and Hawking [2] as the four laws of black hole mechanics, in analogy with ordinary 
thermodynamics. They showed that (i) the surface gravity k, which is the force applied by an observer at 
spatial infinity to hold a particle of unit mass in place at the location of the horizon, is same everywhere 
on the horizon for a stationary black hole — the statement of the zeroth law of black hole physics, (ii) The 
surface gravity k of the black hole analogically resembles the temperature (7h) of the hole, in accordance 
with the interpretation of the horizon area as the black hole entropy. This may be conceived from the 
first law, which states that the change in mass (energy) of the black hole is proportional to the surface 
gravity times the change in horizon area. 



* email: saurya.das@uleth.ca 

t email: shanki.subramaniam@port.ac.uk 

* email: sourav. sur@uleth.ca 



1 



Hawking's demonstration of black hole thermal radiation [3] paved the way to understand the physical 
significance of the temperature Th (and hence the entropy-area proportionality). Hawking showed that 
quantum effects in the background of a body collapsing to a Schwarzschild black hole leads to the emission 
of a thermal radiation at a characteristic temperature: 

hc\ k ( h<? \ 1 



Th \k B ) 2ir \Gk B )^M , (1) 

where G is the Newton's constant in four dimensions, k B is the Boltzmann constant, and M is the mass 
of the black hole. The factor of proportionality between temperature and surface gravity (and as such 
between entropy and area) gets fixed in Hawking's derivation [3] , thus leading to the Bekenstein-Hawking 
area law: 

'k B \ An 



s m - ( -f ) ~pT ' ( 2 ) 



where £ pl = \J Gfij(? is the four dimensional Planck length. 

Black-hole thermodynamics and, in particular, black-hole entropy raises several important questions 
which can be broadly classified into two categories: 

• Gravitational collapse leading to black-hole formation 

(i) What is the dynamical mechanism that makes S BH a universal function, independent of the 
black-holes past history and detailed internal condition? 

(ii) How does a pure state evolve into a mixed (thermal) state? Is there a information loss due to 
the formation of black-hole and Hawking process? Does the usual quantum mechanics need to 
be modified in the context of black-holes? 

(iii) Can quantum theory of gravity remove the formation of space-time singularity due to the 
gravitational collapse? 

• Near thermodynamical equilibrium 

(iv) Unlike other thermodynamical systems, why is black-hole entropy non-extensive? i. e. why 
S BH is proportional to area and not volume? 

(v) Why is the black-hole entropy large? 1 

(vi) How S BH concords with the standard view of the statistical origin? What are the black-hole 
microstates? 

S = k B In (# of microstates) 

(vii) Are there corrections to S BH ? If there are, how generic are they? 

(viii) Where are the microscopic degrees of freedom responsible for black-hole entropy located? 



1 In order to see that, let us compare S Bli with the entropy of the current universe. The entropy of the universe within 
our horizon today [4] is 



\2.75K J 



5 Unlv ~10 87 (^— ) (3) 



where T is the temperature of the universe, h is of the order unity. On the other hand the entropy of a Schwarzschild 
black-hole is 

s..~>°"G£) 

where Mq is the solar mass. Hence, a couple of hundred thousand solar mass black holes can contain as much entropy as is 
free in the entire universe. There is increasing evidence that super-massive black holes exist at the center of many galaxies. 
By now we know that a large fraction of galaxies — of the total (of order) 10 11 — contain such super-massive black holes 
with mass range 10 6 < M bh /Mq < 10 10 . This implies that the entropy of the black holes dominates all other sources of 
entropy. Hence, understanding the origins of black-hole entropy may help us explain the entropy budget of the universe. 
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These questions often seem related, which a correct theory of quantum gravity is expected to address. 
In the absence of a workable theory of quantum gravity, there have been several approaches which address 
one or several of the above questions. Most of the effort in the literature, as in this review, has been 
to understand the microscopic statistical mechanical origin of S BH assuming that the black-hole is in a 
(near) thermal equilibrium or not interacting with surroundings. 

The various approaches may broadly be classified into two categories [5]: (a) the ones that associate 
S BH with fundamental states such as strings, D-Branes, spin-networks, etc. [6,7], and (b) the other that 
associate S BH with quantum fields in a fixed BH background, like the brick- wall model [8], the quantum 
entanglement of modes inside and outside of the horizon [9-12] and the Noether charge [13]. As mentioned 
above, although, none of these approaches can be considered to be complete; all of them — within their 
domains of applicability — by counting certain microscopic states yield (2). This is in complete contrast 
to other physical systems, such as ideal gas, where quantum degrees of freedom (DOF) are uniquely 
identified and lead to the classical thermodynamic entropy. 

The above discussion raises three important questions which we try to address in this review: 

1. Is it sufficient for an approach to reproduce (2) or need to go beyond S BH ? 

As we know, S BH is a semi-classical result and there are strong indications that Eq. (2) is valid 
for large black holes [i.e. Ah 3> ^ P1 ]- However, it is not clear, whether this relation will continue 
to hold for the Planck-size black-holes. Besides, there is no reason to expect that S* BH to be the 
whole answer for a correct theory of quantum gravity. In order to have a better understanding of 
black-hole entropy, it is imperative for any approach to go beyond 5 BH and identify the subleading 
corrections. 

2. Are the quantum DOF that contribute to S Bli and its subleading corrections, identical or different? 

In general, the quantum DOF can be different. However, several approaches in the literature [14] 
that do lead to subleading corrections either assume that the quantum DOF are identical or do not 
disentangle DOF contribution to 5 BH and the subleading corrections. 

3. Can one locate the quantum DOF that give rise to S BH and its subleading corrections? 

More specifically, can we determine to what extend do the quantum DOF close to the horizon or far 
from the horizon contribute to the S BH and its corrections. Depending on the approach, one either 
counts certain DOF on the horizon, or abstract DOF related to the black hole, and there does not 
appear to be a consensus about which DOF are relevant or about their precise location [5]. 

In this review, using the approach of entanglement of modes across the black-hole horizon, we address 
the above three issues and show that: 

1. the entanglement leads to generic power-law corrections to the Bekenstein-Hawking entropy, of the 
form [15]: 

s '= ao (i?) + ai (i?) " ; ( 0<i,<1 )' ( 5 ) 

where so,a\ are constants and a is the lattice spacing, 

2. the quantum degrees of freedom that lead to S BH and subleading corrections are different, and 

3. the contribution to S BH comes from the region close to the horizon while the subleading corrections 
have a larger contribution from the region far from the horizon [16]. 

This review is organized as follows: in the next section we briefly review the basic features of quantum 
entanglement and the concept of entanglement entropy. In sec. (3), we provide a heuristic picture of the 
link between the entanglement entropy and black-hole entropy. In sec. 4 we discuss the procedure and 
assumptions to compute the entanglement entropy of a scalar field in black-hole space-times. In sec. 5 
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we review the cases where the scalar field is either in its ground state, or in a generalized coherent state, 
or in a class of squeezed states — for all these cases the area law is found to hold. In sec. 6, we show 
that the for the superposition of ground and 1-particle state, the entanglement entropy has a subleading 
power-law corrections to the AL. In sec. 7 we study the locations scalar field DOF that are responsible for 
the entanglement entropy. In sec. 8 we examine the entanglement entropy due to a massive scalar field 
and compare it with that obtained for a massless scalar field. We conclude with a summary and open 
questions in sec. 9. In Appendix A we discuss the relevance for considering massless or massive scalar field 
for computing the entanglement entropy of black holes, from the perspective of gravitational perturbations 
in static black hole space-times. In Appendix B we discuss the steps to obtain the Hamiltonian of a scalar 
field in the static black-hole background, which in Lemaitre coordinates, and at a fixed Lemaitre time, 
reduces to the scalar field Hamiltonian in flat space-time. 

Before we proceed, we outline our conventions and notations: We work in four-dimensions and our 
signature for the metric is (— , +, +, +). Hereafter, we use units with k B = c = h = 1 and set the Planck 
mass Mp = 1/(16^(7). The quantum field ip is a minimally coupled scalar field. 

2 Entanglement entropy 

Let us consider a bipartite quantum mechanical system, i.e., a system which can be decomposed into 
two subsystems u and v, such that the Hilbert space of the system is a tensor product of the subsystem 
Hilbert spaces: 

n = n u ®n v . (6) 

Let \ui) and \vj) are eigen-bases which span the Hilbert spaces TL U and TL V respectively. Then \ui) ® \vj) 
forms an eigen-basis in TC, and in terms of this a generic wave-function in TL can be expanded as: 

|*> = XX-|«i> ® !«,■> £7i - w 

ij 

The density matrix operator of the whole system defined by: 

P=l*><*|, (8) 

has the following properties: (i) it is non-negative, i.e., for any vector \cf>), {(p\p\4>) > 0, (ii) it is self-adjoint 
(pt = p), (iii) its trace is unity (Tr(p) = 1), and (iv) it is idempotent (p 2 = p), which means that it has 
only two eigenvalues p n = 0, 1. 

Now the expectation values of the operators which represent the observables in one subsystem, say u, 
can be obtained using the reduced density matrix operator (p u ) for the subsystem u, which is the trace of 
the full density matrix operator (p) over the other subsystem v: 

p u = Tr v (p) = J2(vi\p\vi) = Wijdik^u^ . (9) 

l i,j,k 

When the correlation coefficients dij = 1 (for all the wave-function describing the entire 
system is just the product of the wave-functions \^ u ) and \^ v ) which describe the subsystems u and v 
respectively, and the full density matrix is p = p u p v . The reduced density matrix for the subsystem u 
(say) is then simply given by p u = \^ u ){^ u , and bears the same properties as those listed above for the 
full density matrix. In general, however, |*) / (g) and the two subsystems are said to be in 

an entangled or EPR state 2 . For such entanglement the reduced density matrix (p u ) for the subsystem 
u is given by the general expression (9). A similar expression can be obtained for the reduced density 
matrix (p v ) for the subsystem v, by taking the trace of full density matrix (p) over the subsystem u. 

2 For details, see the review [17]. 
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Both p u and p v have the same properties as listed above for p, except that the idempotency is lost, i.e., 
Pu Pu>Pv Pv- As such the eigenvalues of p u and p v are no longer restricted to and 1, but are in 
between and 1. However, an important property of reduced density matrices is that irrespective of the 
tracing over the eigenvalues remain the same, i.e., both p u and p v have the same set of eigenvalues. 

The entanglement entropy or Von Neumann entropy is a manifestation of one's ignorance resulting 
from tracing over one part of the system, and is defined by 

S cnt = -Tr u (p u lnp u ) = -Tr v (p v lnp v ) = -^pnlnpn ; (0 < p n < 1) . (10) 

n 

Since the eigenvalues p n of p u and p v are the same, 5 en t is the same regardless of which part of the system 
is being traced over. 

3 Entanglement entropy and black-hole entropy — Connection 




I — ► r 

Figure 1: Space-time diagram depicting the collapse 
of a star to form a black-hole. 



Before we proceed with the setup and assumption, 
in this section, we provide a heuristic link between 
the entanglement entropy and black-hole entropy. 

To understand this, let us consider a scalar field 
on a background of a collapsing star. Before the 
collapse, an outside observer, atleast theoretically, 
has all the information about the collapsing star. 
Hence, the entanglement entropy is zero. During 
the collapse and once the horizon forms, S BH is 
non-zero. The outside observers at spatial infinity 
do not have the information about the quantum 
degrees of freedom inside the horizon. Thus, the 
entanglement entropy is non-zero. In other words, 
both the entropies are associated with the existence 
of horizon 3 . 

Beside this, 5 ent and 5 BH are pure quantum 
effects with no classical analogue. Based on this 
picture, we obtain S cnt of scalar fields in a fixed 
background. Although, this analysis is semiclassi- 
cal, since the entanglement is a quantum effect and 
should be present in any theory of quantum grav- 
ity, the results presented here do have implications 
beyond the semiclassical regime. 



4 Entanglement entropy of scalar fields — Assumptions and setup 

In this section, we consider a free scalar field propagating in spherically symmetric black-hole space-times. 
The motivation for such scalar field is presented in Appendix A] 

We consider a massless scalar field (cp) propagating in an asymptotically flat, four-dimensional black- 
hole background given by the Lemaitre line-element 93. The Hamiltonian of the scalar field propagating 

3 It is important to understand that it is possible to obtain a non- vanishing entanglement entropy for a scalar field in flat 
space-time by artificially creating a horizon [10]. However, in the case of black-hole, the event horizon is a physical boundary 
beyond which the observers do not have access to information. 
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in the above line-element is 



lm 



oo 



=nf + 

Lm 



(11) 



(12) 



_rVi-/M" im ' V / WW 

where <^ m is the spherical decomposed field and Il; m is the canonical conjugate of (fi m i. e., 

<Plm{r) =r J dn Z lm (9,4>)ip(r) ; Il Zm (r) = r j dVL Z lm {Q , (f>)U(r) , 

where Zi m (9,4>) are the real spherical harmonics. [For discussion about the Lemaitre coordinates and 
explicit derivation of Eq. (11), see Appendix (B).] Note that the above Hamiltonian is explicitly time- 
dependent. 

Having obtained the Hamiltonian, the next step is quantization. We use Schrodinger representation 
since it provides a simple and intuitive description of vacuum states for time-dependent Hamiltonian. 
Formally, we take the basis vector of the state vector space to be the eigenstate of the field operator 
¥>( r >0 on a fixed r hypersurface, with eigenvalues tp(£) i. e. 

^(t,£)M£),t) = </?(OMO^> (13) 

The quantum states are explicit functions of time and are represented by wave functionals \I/[<£>(£), r] 
which satisfy the functional Schrodinger equation: 



i— = J dtH BH (T)*[ip(0,T}. 



(14) 



We, now, assume that the above Hamiltonian evolves adiabatically. Technically, this implies that 
the evolution of the late-time modes leading to Hawking particles are negligible. In the Schroedinger 
formulation, the above assumption translates to \l/[<^(£), r] being independent of time. At a fixed Lemaitre 
time, Hamiltonian (11) reduces to the following flat space-time Hamiltonian [see Appendix (B) for details] 

... . 2 + ^»|. W> 

lm lm 

The quantum states [defined above in Eq. (13)] of this Hamiltonian is time-independent and \I/[<£>] satisfies 
the time independent Scrodinger equation 



I'm J<m I 



JO 



drHV[<p(r)] = 



(16) 



The procedure of finding the entropy involves the following steps: 



(i) Discretize Hamiltonian (15): Obtaining an analytic expression for the entropy using the Von Neumann 
definition (10) is prohibitively difficult. For the field theory, even if we obtain closed-form expression of the 
density matrix, it is not possible to analytically evaluate the entanglement entropy. Hence, we discretize 
the Hamiltonian in a spherical lattice of spacing a such that r — ► rf, r^+i — ri = a. The ultraviolet cutoff 
is therefore M = a~ l .. The lattice is of very large but finite size L = (N + l)a (N 3> 1), with a chosen 
closed spherical region of radius R(n + l/2)a inside of it, as shown in the Fig. 2. It is this closed region, 
by tracing over the inside or outside of which one can obtain the reduced density matrix. We demand 
that the field variables (fi m (r) = for r > L so that the infrared cutoff is M = L _1 .The ultraviolet cutoff 
is therefore M = a -1 . 

Now, in discretizing the terms containing the derivatives in the Hamiltonian (15), one usually adopts 
the middle-point prescription, i.e., the derivative of the form f(x)d x [g(x)] is replaced by fj+i/2[9j+i~9j]/ a - 
The discretized Hamiltonian is given by 



lm 



1 

2~a 



N 

E 

j'=i 



^ImJ + 



3 + 



1 



<Pl 



m,j 



m,. 



•2+1) 



3 + 1 J 



+ 



1(1 + 1) 



<Plm,j 



H 



Y, H > 



lm 



(17) 



lm 
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L=N*a 
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where ipi m j = <pi m (rj), iri m = n m ,j(fj), which satisfy the 
canonical commutation relations: 



[ l Plm,j, 7r l'm',j'] — iSu/S mm i8jji . 



(18) 



Up to an overall factor of o _1 , the Hamiltonian Hi m , given 
by Eq. (17), represents the Hamiltonian of N— coupled 
harmonic oscillators (HOs): 

1 N 1 N 

H (n-ko) = + 2^ x i K V x ii ( 19 ) 

i=l *J=1 

where the coordinates Xi replace the field variables (pi m , 
Figure 2: Discretization of the scalar field the momenta « re P lace the conjugate momentum vari- 
propagating in the flat space-time. The trac- ables and the N x N matnx K *i (hJ = l,-..,N) 

represents the potential energy and interaction between 

the oscillators: 



ing is done over the shaded region. 



K, 



1 



9 / 1 

1(1 + 1) <% + - <5;i<5ji+ 



i + 



1 



+ i 



1 



8t 



ij(i^l,N) 



j(j + 1) 



i + 



l\2 



i(i + 1) 



(20) 



The last two terms which originate from the derivative term in Eq. (15) denote the nearest-neighbor 
interactions. Schematically, the matrix K is as shown below, where the off-diagonal terms represent 
interactions: 



K 



/ x x 

XXX 

XXX 

XXX 

XXX 

XXX 
X X 



(21) 



(ii) Choose the quantum state of the field: The most general eigen-state of the Hamiltonian (19) for the 
iV— coupled HOs is given by 



ip(xi, 



,x N ) 



i=l " ^ 



-k 1/2 x 2 , 



(22) 



where A/i's are the normalization constants given by 



k 1/4 

K Di 



7T 



V 4 ^/2 u n^. 



(i = l,...N), (23) 
,Xjv)j Kd = UKU T is a diagonal matrix with 



x = Ux, (U T U = I N ), x T = (xi, . . .,x N ), x T = (x l5 , 
elements km, and (i = 1 . . . N) are the indices of the Hermite polynomials (7i v ). The frequencies are 
ordered such that koi > koj for i > j. 

The difficulty however is that it is not possible to work with an arbitrary iV— particle state as the 
density matrix (24) cannot be expressed in a closed form. Therefore one has to make some specific choices 
for the quantum state in order to make the calculations tractable. In the following two sections we will 
discuss several possible choices. 
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(iii) Obtain the reduced density matrix by tracing over certain closed region: The reduced density matrix 
is obtained by tracing over the first n of the N oscillators: 



/n 
Y[ dxi Y>(xi, ...,x n ;t) tp*(x 1 , ...,x n ; t') 
i=i 

N 

n^^(^)ex P 



J^J dx 

i=i 



x T £lx 



x fT flx' 



N 

n^(*«*0 ( 24 ) 



J"=l 



where = U T K^ 2 U is an iV x N matrix, such that \Q\ = \Kn\ 1 / 2 , and we have made the following 
change in notation: x T = [x\, . . . , x n ; t±, . . . , ijv-i) = {x±, ■ ■ ■ , x n ; t), with t = t\, . . . ,tN- n ! = ^n+j ; j = 
1, . . . , (N — n). One may verify that p 2 / p, i.e., the state obtained by integrating over n of the HOs is 
mixed, although the full state is pure. 



(iv) Compute the entropy: The entanglement entropy can be calculated by substituting the reduced 
density matrix (24) into the expression (10). 



5 Warm up — Entanglement entropy for (displaced) ground state 

In this section, we review the procedure for obtaining 5 ent in the ground state, coherent state (which is 
a displaced ground state) and a class of squeezed state (which are unitarily related to the ground state). 
We show that the in these three cases that the entanglement entropy is proportional to area. 

Ground State (GS): 

When all the HOs are in their GS, then by setting Vi = 0, for all i, the wave function (22) takes the 
form 



V> gs (xi,...,Xtv) = j ] 



i=l 



7T 



exp 



-A: 1/2 x 2 



— ex P 



7T 



x T ■ Q ■ x 



(25) 



The corresponding density matrix can be evaluated exactly as [10]: 

Pas (*;0 : 
where we have decomposed 



7T 



\n\ 

N-n\ 



exp 



t T ~/t + t /T 7t' 



+ t T f3t' 



(26) 



A B 
B T C 



and defined 



P = 



B T A- 1 B 



c-p. 



(27) 



(28) 



A is an n x n symmetric matrix, B is an n x (N — n) matrix, and C, (3, 7 are all (N — n) x {N — n) 
symmetric matrices. The matrices B and f3 are non-zero only when the HOs are interacting. 
Performing a series of unitary transformations: 



VjV 1 = 1D = diag , (3 = lD 1/2 VpV r lD 



rT-1/2 



WPW T = p D = diag , v = W T ~ffv 



(29) 
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one can reduce p GS (t; t') to a product of the reduced density matrices p (2 _ HO) (t; t') for (N — n) two coupled 
HOs with one oscillator traced over (i.e., N = 2,n = 1) [10]: 



Pas (*;*') 



iV-n 
i=l 



(2-HO) 



(2-HO) 



(t;f) 



|fi| 



exp 



V? + V f 



(30) 



where Vi £ v and /3j G /3. Correspondingly, the total entanglement entropy is a sum of (N — n) two-HO 



entropies S, 



(2-HO) 



N 



S 



(2-HO) 



n) which are obtained using the Von Neumann relation (10) as [10]: 
& , , , Pi 



1-6 



In ^ 



6 



i + a/i-A 



The total GS entropy for the full Hamiltonian H = H-im, Eq.(17), is therefore given by 



(31) 



S GS (n,N) =J2Slm(n,N) = £(2Z + l)5/(n,JV) , Si{n, N) = - ln[l - 



6 



Zm 



z=o 



1-6 



In 6 



(32) 



where (2Z+1) is the degeneracy factor found by summing over all values of m for the Hamiltonian. Ideally 
the upper bound l ma x of the sum should be infinity, however in practice one assigns a very large value of 
Imax to compute the entropy upto a certain precision. The precision limit Pr is set by demanding l max to 
be such that the percentage change in entropy for a change in I by a step size l s t, never exceeds Pr, i.e., 

S(l r nax) S(l max 1st) 



S(l, 



1st) 



x 100 < Pr . 



(33) 



Thus Pr = 0.01 (say) implies that the numerical error in the computation of the total entropy is less than 
0.01%. 

It has been shown in ref. [10] that for N ~ 60 and 
n < N/2, the GS entropy computed numerically 
follows the relation: 



In S„ c vs In (R/a) 



S, 



0.3 ( n+ - 



:r 

0.3 - 



0.3 

47T 



A 



(34) 



In S = 1.9904 In (R/a) - 1.2097 



where R = (n + l/2)a is the radius of the hypo- 
thetical spherical surface (which is an analogue of 
the horizon in flat space-time) the DOF inside or 
outside of which are traced over, and A = AttR 2 
is the surface area of this sphere. For a fixed 
value of n, if N is varied the GS entropy re- 
mains the same. This implies that S GS does not 
depend on the infrared cutoff M = L~ l , where 
L = (N + l)a, although it depends on the ultravi- 
olet cutoff M = a~h S GS = 0.3M 2 R 2 [10]. 

In Fig. 3 we have plotted ln5 GS vs R/a = 
n + 1/2 for the parametric values: N = 300 and 
n = 100 - 200. The data fits into a straight line Fi S ure 3: Plot of iogarithm of the GS entropy versus 

, i o /r ■ • v -i. ln(R/a), where R = (n+ l/2)a radius of the hypothetical 

with slope very close to 2 (tor a precision limit , ' ',' v „ \ T „„„ „„„ „„„ ™ 




4.9 5 

In (R/a) 



Pr = 0.1% in the computation performed using 
MATLAB has been). The area law (S GS oc A) is 
thus found to hold generically for GS. 



sphere (horizon), for N = 300, n = 100 - 200. The graph 
is a straight line with slope ~ 2. 
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Generalized Coherent State (GCS): 



We now examine the situation where all the oscillators are in GCS, instead of being in GS. Similar to 
GS, a GCS is a minimum uncertainty state {AxAp = 1/2). However unlike GS, the GCS is not an energy 
eigenstate of HO, it is an eigenstate of the HO annihilation operator with real eigenvalue. The GCS 
wave- function V'gcs differs from the GS wave-function, Eq. (25), by constant shifts oti in the coordinates 
Xj, as: 



ip GCS (xi,...,x N ) 



N , , \ 1/4 

n v exp 



N 



i=l 



1 



F[ exp (-ipjOi) exp ( -^fc^ 2 xf 



(35) 



where Pj = —id/d^ are the momenta. Physically, the real and imaginary parts of the N complex 
GCS parameters oti correspond to the classical position (xo) and momentum (po) of the individual HOs 
respectively, i.e., ctij = xq — ipo/hpi- Defining shifted coordinate variables [11]: 



x = x — U 1 a 



dx = dx . 



one can show that 



ip GCS (xi,...,x N ) 



/|0| N 


V4 


x T • Q • x 




) exp 


2 



ip GS (xi,...,x N ) 



(36) 



(37) 



Thus the GCS wave-function is of the same form as the GS wave-function, albeit in terms of the shifted 
variables. As such, the reduced density matrix p GCS for GCS will also be of the same form as that for GS, 
Eq.(26), with the variables (t; t') replaced by (t; t'), where t = t±, . . . ,tN-n ', tj = x n+ j, j = 1, . . . , (N — n): 



/n 
Y[dxi V> GC s( x i' ■■■,x n ]t) ip GCS (xi, . . .,x n ;t') = p GS (t;t') . 
i=i 



(38) 



Consequently, the eigenvalues of the density matrix, and hence the entropy for GCS will be the same as 
those for GS [11]. Thus the area law holds for the case of GCS as well. 

Squeezed State (SS): 

Let us now consider the case where all the HOs are in a class of SS. The squeezed states are also 
minimum uncertainty packets (AxAp = 1/2), similar to GS and GCS. However, there is a squeezing 
either in the positions (Ax <C 1, Ap » 1) or in the momenta (Ap <C 1, Ax » 1). We consider a class of 
SS, which is characterized by a unique squeezing parameter £, and described by the wave-function 



v 



ii> ss (xi,...,x N ) ~ c^n^p 



i=i 



N 



1=1 



In this case, defining scaled coordinate variables [11]: 

x = y^C x , dx = -y/C dx , 



(39) 



(40) 



one finds that the SS wave-function and the corresponding density matrix reduce to the forms same as 
those for GS, albeit with the replacement x — > x. As such, upto an irrelevant multiplicative factor, the 
SS entropy turns out to be the same as the GS entropy, Eq. (34), and hence the area law holds. 

Thus, in all the above cases where we have the states which form minimum uncertainty packets, the 
entanglement entropy is proportional to the area. In the next section, we show that if we consider a higher 
excited state, the area law is not robust. In particular, we show that the superposition of the ground and 
excited state leads to power-law corrections to area. 
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6 Power-law corrections to the area-law 



Ideally the most interesting thing would be to find the entanglement entropy for a general eigenstate, given 
by Eq. (22), which may be the ground state or the first, second, etc. excited states, or a superposition 
of such states. However, as mentioned earlier, the density matrix cannot be written in a closed form for 
such a general state. We therefore resort to a simple class of excited states, viz., (i) 1-particle excited 
state and (ii) superposition of GS and 1-particle excited state. In the case of 1-particle excited state, we 
show that the entanglement entropy does not scale as area while in the superposed state the entanglement 
entropy has power-law corrections to the area-law. 




6.1 1-Particle Excited state 

1-particle excited state is described by a wave-function ^ ES which is a linear superposition of N HO wave 
functions, each of which has exactly one HO in the first ES and the rest (N — 1) in their GS [11,12]. Using 
Eq. (22), such an 1-particle ES wave-function is expressed as: 

^ ES (xi,...,x N ) = Yl (~^) aiHl { k Di^i) ex P (^~\ Yl k Dj 

= V2(a T K 1 J\^ GS (x 1 ,...,x N ), (41) 

where a T = (a>i, . . . , ajv) are the expansion coefficients, and a T a = 1 so that i[) ES is normalized. 
From Eq.(24) one finds the reduced density matrix p ES (t; t') for ES as: 

/n „ n 

\{ dx^ ES (xi- 1) V ES fa; t') = 2 H d Xl [x' T Ax] V GS (x i; t) V GS {xf, t') , (42) 
i=i J i=i 

where A is a N x N matrix given by 



A^. If K}{' a a? & U = ( ^ ^ ). 



A^4 is an n x n symmetric matrix; Ag is an n x (N — n) matrix; Ac is an (N — n) x {N — n) symmetric 
matrix. Defining two {N — n) x (N — n) square matrices A^ and A 7 such that 

A/3 = ^(2A C - A T B A- l B - B T A- 1 A B + B T A- 1 A A A~ 1 B) 

A 7 = i(2A^- 1 5 - B T A- 1 A A A~ 1 B) , (44) 

K 

one can express the above density matrix for ES in terms of the density matrix for GS as [11]: 

1 



Pes (*;0 = 



1 _ t T A 1 t + t' T A 1 t' +t T A/3t> 



Pas (*;0, (45) 



where k = Tt^aA^ 1 ). In the above, the matrix A^ is symmetric, whereas the matrix A 7 is not necessarily 
symmetric due to the presence of the first term in the parentheses on the right hand side. 

Eq. (45) is an exact expression for the density matrix for a discretized scalar field with any one HO 
in the first ES and the rest in the GS. However, unlike the GS density matrix p GS , given by Eq. (26), 
the ES density matrix p ES , Eq. (45), contains non-exponential terms and hence can not be written as a 
product of (N — n), 2-coupled HO density matrices. Therefore, the ES entropy cannot be expressed as a 
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sum of 2-HO entropies, as in the case of GS. One may however note that the GS density matrix p GS , Eq. 
(26), is a Gaussian that attenuates virtually to zero beyond its few sigma limits. Therefore, if 



£l — t max Apt max 1 , 62 — ^max-^y^rnax ^ 1 



where 



3(iV - n) 



(1,1, 



max \y/2Tr(r, - 0) 
corresponding to 3<r limits of the Gaussian inside p GS , then one may approximate 

t T A^t + t' T A^t' 



1 



+ t T A t' 



exp 



+ t T A p t' 



(46) 
(47) 

(48) 



Consequently, with a shift of parameters: /?' = + A^, 7' = 7 + A 7 , the ES density matrix above can also 
be approximated as a Gaussian 



/ A 1 



(49) 



In S ES vs In (R/a) for = 20 



lnS ES vsln(R/a)foro = 30 



Factorizing this once again into N — n two-HO density matrices, albeit in terms of the shifted parameters 
(/?', T 7 ), one can evaluate the ES entanglement entropy S ES . 

The entropy computation is done nu- 
merically (using MATLAB) in [11], with a 
precision setting Pr = 0.1% for the set of 
parametric values: N = 300, n = 100 — 
200, o = 10 — 50, where o is the number 
of last non-vanishing entries in the vector 
a T , i.e., a T = (1/^/5) (0, •• ■ ,0;1,-- - ,1). 
The criteria (46) are found to be satisfied 
for these choices of the parameters. The 
InS'gg vs. ln(i?/a) data for different fixed 
values of the amount of excitation o (= 
20, 30, 40, 50) fit approximately to straight 
lines as shown in Fig. 4. The slopes of 
all these straight lines are less than 2 and 
the higher the value of o the smaller is the 
slope. This shows that the ES entropy ap- 
proximately scales as a power of the area: 




4.7 4.8 4.9 5 5.1 
In (R/a) 

In S ES vs In (R/a) for = 40 



In S ES = 1.7695 In (R/a) - 0.0337 



4.7 4.8 4.9 5 5.1 5.2 

In (R/a) 

lnS ES vsln(R/a)foro = 50 



Ao 
4tt 



; A 



AttR 2 . 
(50) 




In S ES = 1.6510 ln(R/a) + 0.6158 



8.5 




4.7 4.8 4.9 5 

In (R/a) 



5.1 5.2 



In S ES = 1.5218 ln(R/a) + 1.3498 



4.7 4.8 4.9 5 5.1 5.2 

In (R/a) 



Figure 4: Plot of logarithm of the ES entropy versus ln(i?/a) for 
N = 300, n = 100 - 200 for different values of o = 20,30,40,50. 
The data approximately fit to straight lines with slope < 2. 



The power //, however, is always less than 
unity and decreases with the increase in 
the number of excitations o. The coeffi- 
cient Ao on the other hand increases with 
o. Thus contrary to the cases of GS, GCS 
and SS, the AL is always violated in the 
case of the 1-particle ES [11]. 

Now, it is quite interesting to see that even for a small amount of excitation (o ~ 20) the entropy-area 
relationship is so drastically changed that the AL could not be recovered in any limit. Further studies 
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with computations of higher precision (Pr = 0.01%) [15] however reveal that there are slight variations 
in the linear fits of hiS ES vs. ln(i?/a) data for different values of o. In fact, it has been shown that the 
ES entropy actually approaches the GS entropy (and hence obeys the AL) for very large area [15]. Thus 
the excitations seem to give rise to some corrections to the AL which are significant for smaller areas, but 
become negligible for very large areas. These corrections are not manifested by the linear fits of the data 
(shown in Fig. 4). Therefore a more accurate non- linear fitting is required. We will discuss about these 
corrections to the AL in the next section where we consider the scalar field to be in a more general state 
— the MS — which is a linear superposition of GS and 1-particle ES. 

6.2 Superposition of ground and excited state 

The superposition of ground and excited state wave-function ip MS is given by 

V> MS 0M) = [co^ G sOM) + ci V'es (£;*)] . (5i) 

where if) GS * s the GS wave-function, given by Eq. (25), V E s ^ s the ES wave-function, given by Eq. (41), 
x = {xi, ■ ■ ■ ,x n } , and as before tj = x n+j (j = 1, ■ ■ ■ , N - n) ; t = {h, ■ ■ ■ ,t N ^ n } = {x n+1 , ■ ■ ■ ,x N }. 
We assume that cq and c\ are real constants, and ^ MS is normalized so that Cq + c\ = 1. Using Eq. (41), 
we can write, 

V'msOM) = [co + ci f(x;t)]i/j GS (x;t) , f(x;t) = \/2a T K^Ux = y T x, (52) 

where the column vector a includes the expansion coefficients defined in the previous section [a T = 
(«i, . . . , aw) = (l/ v / o)(0, . . . , 0; 1, . . . , 1)], and y is an iV-dimensional column vector y defined as 

y = V2U-K^a = ( £ ) (53) 

yA and ys are n- and (N — n)-dimensional column vectors, respectively. 
The MS density matrix can be expressed as a sum of three terms: 

/n 
Y[dxi ^ MS (x;^* s (£;*') = c 2 p GS (t;t>) + c\ p ES (t;t') + Cod p x (t; t') , (54) 
i=i 

where p GS (t;t') is the GS density matrix, Eq. (26), p ES (t;t') is the ES density matrix, Eq. (42), and 
p x (t; t') is a cross term due to the superposition of GS and ES. Identifying the matrix A, its components, 
and the constant k (defined for ES in the previous section), with the column vector y and its components: 



i T 
-jVBVB , 

(55) 



(VB-pf (t + f) p GS (t;t'), (56) 



(57) 



A= 2 Vy = US A c) ' 

= ^vava ; a b = \vavb ; A c 

k = Tt{K a A- 1 ) = X -y T A A- x y A , 
one can evaluate the cross term p x as 

„ n 

Px (t; t') = / J] dx i l ) + V'gs *) V> 



1=1 



where p is an (N — n)-dimensional column vector defined by 

p = B T A- 1 y A - 
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The full density matrix for the MS, Eq. (54), thus reduces to 

PM S (t;t')= [cl + clK{l + U (t;t')}+c c 1 v(t;t f )] p GS (t;t') , (58) 



where 



Defining 



where 



fT A f _i_ flT A f I 

= _ A^t + t Ajt + ^ ^ = (yB _ p) T (t + ^ (59) 



K 2 

F(t;0 = 1 + Kiwftt') + K2«(t;0 + -^-v 2 (t;t'), (60) 



w(t;t') = J — ?— + fkpd', 

A/3, = A/3 - 2k (Ap - , Ay = A 7 + 2k ^ - ^ 

C C l C C 1 -2,2 

K = — 5 Kl = — ; K2 = —=~ , K = Cn + CiK, (61) 

K K K 

the density matrix (58) can be written as 

p MS (t;t') = RF{t-t>) P GS (*;*')■ (62) 

In the above, Ap/ and Ay are (N — n) x (N — n) matrices, and constants (kq, ki, k 2 ) describe the amount 
of mixing between the GS and ES. 

Now, similar to the case of ES, here also the pre-factor F(t;t') of the Gaussian p GS (t;t') in Eq. 
(62) contains non-exponential terms. Therefore p MS cannot be factorized into (N — n) two-HO density 
matrices. However, as mentioned in sec. 6.1, that when the vector t T is outside the maximum i^ aa; , 
Eq. (47), corresponding to the 3cr limits, the Gaussian inside p GS (t;t') is negligible. Therefore, if the 
conditions (46) (given in sec. 6.1) as well as the conditions 

h = tmax^P'tmax < 1 , e 2 = i^ aa; Ay i maa; < 1 (63) 

are satisfied, then keeping terms up to quadratic order in t, t', the pre-factor F(t; t') can be approximated 
as 

F(t;t') » exp [ Kl w(t;t') + k 2 v(t;t')] . (64) 
Then by using Eq.(26) for p GS (t;t') the (approximated) MS density matrix can be expressed in the form: 



PmsM = fty ^J-S^i exp[z(t^) + K2v(t;l/)] , z(t;t') = - ^ ± - ^ 7 - + t T (3't' , (65) 
where the (N — n) x (N — n) matrices f3' and 7' are defined by 

0' = + K\Api = + K\A/3 - 2k K! (a^ - ^pj 

7' = 7 + K]Ay = 7 + KiA 7 + 2k Ki ^A/3 - ^-^j . (66) 

The matrix f3' is symmetric while the matrix 7' is not necessarily symmetric. 

Now shifting (N — n) variables t = {x n +i,--- ,xn} and t' = {x' n+1 ,--- ,x' N } by constant values 
s = {si,-- - ,s N - n }: 

*-►* + «; t' -► t' + s , (67) 
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the density matrix (65) becomes 



p MS (*;0 = M ex P 



t T it + t ,T it' 



+ t T f3't' 



M = R\ 



7T 



JV-n 



\A\ 



cxp 



The (iV — n)-dimensional constant column vector s is determined from the equation 



y 7' + 7' 



(68) 



(69) 



For either Co = or ci = 0, the constant K2 = 0, whence from the above equation (69), we have s = 0. 
Thus it can be verified that the MS density matrix (68) reduces to the GS density matrix (26) when 
Co = l,ci = 0, which imply (3 1 = (3,j' = 7. On the other hand when cq = 0, c\ = 1, which imply 
(3' = f3 + A/3, 7' = 7 + A 7 , the MS density matrix (68) is the same as ES density matrix (42). In general, 
when both cq and c\ are non-vanishing, then under the shifts (3 — > (3', 7 — > 7' (where /?' and 7' are given 
by Eqs. (66)) the MS density matrix (68) is of the same form as the GS density matrix (26), up to a 
normalization factor given above. Such a normalization constant does not affect the entropy computation. 
Therefore the total MS entropy 5 MS can be evaluated following the same steps [Eqs. (29) - (32)] as in 
the case of GS, albeit with the replacements (3 — > (3', 7 — ► 7'. 

Computation of the entanglement entropy has been done numerically (using MATLAB) in [15], for 
N = 300, n = 100 - 200, o = 30, 40, 50, and with a precision setting Pr = 0.01% in each of the following 
cases: 

(i) GS (c = l,ci =0), 

(ii) ES ( Co = 0,ci = 1), 

(iii) an equal mixing (MS^g) of ES with GS (co = c\ = l/y/2), and 

(iv) a high mixing (MS Hi ) of ES with GS (c = 1/2, c x = y/3/2). 

The conditions (46) as well as (63) are found to be satisfied for the above values of the parameters. 

Before proceeding to the results we would like to mention the following: the expectation value of 
energy, 6, for MS turns out to be 



N 



£ = (4; MS \H\i; MS }=£ + j- k 



1/2 

Di ' 



(70) 



i=N-o+l 



where £ q = \ YliLi ^th ls the (zero-point) GS energy. The fractional excess of energy over the zero-point 
energy is therefore given by 



A£ 
£0 



£-£ n 



24 



sr^N-o . 1/ 



1 + 



1/2 

i 



E 



N 

i=N-o+l "'Di 



1/2 



(71) 



Now, the value of c\ is between and 1 and as mentioned earlier k^i > kpj for i > j. Therefore even 
in the extreme situation c\ = 1, i.e., ES, with a fairly high amount of excitation o ~ 50, the fractional 
change in energy is at most about ~ 4%. Moreover, since there are o number of terms in the sum in the 
second term of Eq.(70), the excitation energy (£ — £$) ~ 1 (in units of 1/a, where a is the lattice spacing). 
Hence, if we choose a ~ Planck length, then this excitation energy is of the order of Planck energy. As the 
mass of a semi-classical black hole is much larger than the Planck mass, one may therefore safely neglect 
the back-reaction of the scalar field on the background. 

Fig. 5 shows the plots of the logarithm of the total entropy S vs. ln(i?/a) = ln(n + 1/2), for the cases 
of GS, ES and MS (Eq/Hi) with different values of the excitation (o = 30, 40, 50). The plot for GS is the 
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In S vs In (R/a) plots for GS, ES & MS(Eq/Hi) 




In (R/a) In (R/a) In (R/a) 

Figure 5: Plots of logarithm of GS, ES and MS (Eq/Hi) entropies versus \n(R/a) for N = 300, n = 100 - 200 and 
o = 30, 40, 50. The numerical error in the computation is less than 0.01%. 
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Asymptotic behaviour for o = 50 
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Figure 6: Ratios of GS and MS (Eq/Hi) or ES entropies and their reciprocals plotted against the area A (in units 
of a 2 ) for o = 30,50, to show the asymptotic nature of the MS and ES entropies with respect to the GS entropy. 
The curves on the upper half (above 1) show the variation of S XS /S GS with A, where XS stands for MS (Eq/Hi) or 
ES, while the lower curves show the variation of S GS /S XS with A. 



linear fit with slope ~ 2 shown earlier. The plots for the MS (Eq/Hi) cases, as well as for ES, are nearly 
linear for different values of the excitations o = 30, 40, 50 and appear to coincide with the plot for GS 
for large areas (.4 = 4ttR 2 S> a 2 ). For a closer examination of this, the ratios 5 MS (EqorHi)/S' GS , S ES /S GS 
and their inverse are plotted against the area A in Fig. 6. All these ratios approach to unity with 
increasing area for different excitations (o = 30,50), i.e., the MS (Eq/Hi) and the ES entropies coincide 
asymptotically with the GS entropy, following the criterion of 'asymptotic equivalence' [18]: 



lim ^—t^ = 1 5 lim ^tj^ = 1 ' XS = MS (Eq or Hi) or ES . (72) 



From Fig. 6 one can also observe that the MS(Eq) entropy is closer to the GS entropy for large A, than 
the MS (Hi) entropy and the ES entropy, the latter being the farthest. Thus for smaller values of the 
relative weight c% of the mixing of ES with GS the asymptote is sharper. 

Fig. 7 shows the best fit ratios of the MS entropies (for equal and high mixings, with o = 30, 40, 50) 
to the GS entropy, which follow a simple formula: 

t^' + ^i^T- <73> 
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Figure 7: Best fit plots (solid lines) of the relative mixed state entropies (S MS /S G3 ) for equal and high mixings 
versus the area A (in units of a 2 ), for o = 30, 40, 50. The corresponding data are shown by asterisks. 



Fitting 
Parameters 


For MS Bg 


For MS m 


o = 30 


o = 40 


o = 50 


o = 30 


o = 40 


o = 50 


°o 


1.001 


1.002 


1.003 


1.001 


1.004 


1.006 




1738 


4288 


8039 


2956 


7652 


14120 


V 


1.180 


1.210 


1.225 


1.141 


1.178 


1.192 



Table 1: Values of the parameters of the fit S ms /S GS = ao + cri (-4/a 2 ) v for both MS(Eq) and MS(Hi) cases with 
excitation o = 30, 40, 50. 



where the values of the fitting parameters &q , u\ and v are shown in Table 1 for different values of 

= 30,40,50. For all these values of o, the parameter <5"o ~ 1 in both MS(Eq) and MS(Hi) cases. The 
parameter a± is of the order of 10 3 and increases with increasing excitations. The parameter v lies between 

1 and 1.25 for the above values of o, and also increases with increasing o. Using the expression for the 
GS entropy, viz., 5 GS = no{A/a 2 ), where no is a constant, we can rewrite the above Eq. (73) as 

5 ms = °~o + o-i (J^J , (74) 

where oq = rio<5"o,o"i = n^ai oc c\ and v = v — 1. The exponent —v lies between and —0.25 for both 
equal and high mixings with the above values of o. It is instructive to stress the implications of the above 
result: 

(i) For the pure vacuum wave-functional, c\ = and S ent is identical to Bekenstein-Hawking entropy. This 
clearly shows that the entanglement entropy of ground state leads to the area law and the excited states 
contribute to the power-law corrections. 

(ii) For large black-holes, power-law correction falls off rapidly and we recover S BH . However, for the small 
black- holes, the second term dominates and black-hole entropy is no more proportional to area. Physical 
interpretation of this result is immediately apparent. In the large black-hole (or low-energy) limit, it is 
difficult to excite the modes and hence, the ground state modes contribute significantly to S ent . However, 
in the small black-hole (or high-energy) limit, larger number of field modes can be excited and hence they 
contribute significantly to S ent . 

(iii) The power-law corrections to the Bekenstein-Hawking area law derived here in the context of en- 
tanglement of scalar fields have features similar to those derived in the case of brick-wall model [19] and 
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higher-derivative gravity [13]. For instance, it was shown that the entropy of five-dimensional Boulware- 
Deser black-hole [20] is given by 

S = j + c A 1/3 ; c = constant . (75) 

As in Eq. (74) the above entropy is proportional to area for large horizon radius, however it strongly 
deviates in the small horizon limit. It is important to note that the corrections to the black-hole entropy 
are generic and valid even for black-holes in General relativity without any higher curvature terms 4 . 

7 Location of the degrees of freedom 

In the previous section, we showed that the entanglement entropy provides generic power-law corrections 
to the area. We also showed that the quantum DOF that contribute to the area law and the subleading 
corrections are different. As we had mentioned in the introduction, this leads to another question: To 
what extend does the quantum DOF close to or far from the horizon contribute to the black-hole entropy? 
In this section, we address this question and show that large contribution (~ 97%) to the area-law comes 
from close to the horizon while the subleading contributions has a larger contribution from the regions 
far from the horizon. 

Let us recall the expression for the interaction matrix K, with elements Kij given by Eq. (20), for the 
system of N HOs. The last two terms which signify the nearest-neighbour (NN) interaction between the 
oscillators, are solely responsible for the entanglement entropy, i.e., if these two terms are set to zero the 
entropy vanishes. In order to find which DOF give rise to the entropy or what are their contributions, let 
us perform the following operations on the matrix K [16]: 

Operation 1: 

Let us set the off-diagonal elements of K, which signify the NN interactions, to zero (by hand) 
everywhere except in a 'window', whose center is at a point q. The indices i,j of the matrix K run from 
q — s to q + s, where s < q, so that the interaction region is restricted to a width of d = 2s + 1 radial 
lattice points. For instance, with s = 1 the window is of size 3x3, and the matrix K is schematically 
depicted as: 



/ x 



K 



|x 

x 



X 
X 
X 



(76) 



x I 
X I 



V 



Now we let the center q of window to vary from to a value greater than n, and thus allow the window 
to move rigidly across from the origin to a point outside the horizon at n as shown pictorially in Fig. 8. 
For every window location (values of q) and fixed window width d, we compute the entanglement entropy 
S(q, fixed d) in each of the cases GS, ES and MS (Eq/Hi) and find the percentage contribution of the 
entropy as a function of q: 

S(q, fixed d) 



pc(q) 



Stotal 



X 100, 



(77) 



4 In this context, it should be mentioned that it is not possible to check for logarithmic corrections to the entropy in our 
analysis, as the numerical error we obtain is much larger than ln(n + 1/2). 
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where Stotai is the total entropy with all NN interactions, i.e., with the indices i,j running from to N. 
The variations of pc(q) with q for a window width of d = 5 lattice points is shown in Fig. 8, for fixed 
values N = 300, n = 100 in each of the cases GS and ES, MS (Eq/Hi) with o = 30, 50. 



GS torn = 100 ESforn = 100 




Figure 8: The figure on the left is a pictorial representation of the entropy due to changing positions of the 
interaction window (green dotted boxes). Contributions to the entropy are observed only when the window includes 
the horizon at n. The figures on the right are the actual plots of the percentage contribution pc(q) to the total 
entropy as a function of window position q, for a fixed window size d = 5 and fixed N = 300, n — 100, in each of 
cases of GS, ES and MS (Eq/Hi). For ES and MS (Eq/Hi) the solid curve is for o = 30 whereas the broken curve 
is for o = 50. 

In all the cases of GS, ES and MS (Eq/Hi) there is no contribution to the entropy, i.e., pc(q) = 0, 
when the interaction window does not include the horizon at n. The contributions to the entropy rise 
significantly when the window includes the horizon, i.e., when the window center q is very close to n. In 
the case of GS, pc(q) peaks when the window is symmetrically placed between inside and outside, i.e., 
when q = n, and decreases when q moves away from n. In the ES and MS(Eq/Hi) cases, however, the 
peak tends to shift towards a value q > n and the amplitude of the peak diminish with increasing o and/or 
the mixing weight c\. As a result, the entire profile of pc(q) is more and more asymmetric with increasing 
o and/or ci, when the window is symmetrically placed between inside and outside of the horizon. The 
peak is shortest and the profile is most asymmetric for ES with o = 50, as shown in Fig. 8. On the whole 
the above results thus confirm that the entanglement between the scalar field DOF inside and outside the 
horizon gives rise to the entropy, and the DOF in the vicinity (inside or outside) of the horizon contribute 
most to the total entropy while the DOF that are far from the horizon contribute a small portion that 
remains. Such contributions from the far-away DOF increase with increasing excitations and/or amount 
of mixing of ES with GS. This is indicated by the diminishing maxima of pc(q) in the MS (Eq/Hi) and 
ES cases, compared to the case of GS. To estimate the contribution of the far-away DOF on the total 
entropy we perform an alternate operation on the interaction matrix K as described below. 

Operation 2: 

let us again set the off-diagonal elements of the matrix K to zero (by hand) everywhere except in a 
'window' whose outer boundary is the horizon at n. The index i of the elements of the matrix K therefore 
runs from a point p to n, where < p < n. Now we let the width d = n — p of the window to vary from 
to n, as shown pictorially in Fig. 9. Computing the entanglement entropy S(d) for every window width 
d, one then finds the percentage contribution of the entropy as a function of d: 

pc{d) = J^- x 100 , (78) 
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where Stotai ^ s the total entropy which is recovered for the full width d = n, i.e., p = 0. For convenience 
we consider here the two extreme cases: GS (o = 0) and ES (with o = 30, 50). The effects of the DOF on 
the entropy in the MS (Eq/Hi) cases are intermediate between those for the GS and ES cases, and are 
therefore of not much interest. 



Entropy 




Figure 9: The figure on the left pictorially shows how the contribution to the total entropy increases with increasing 
width d of the interaction window (green dashed box) in the cases of GS and ES. The figures on the right are the 
actual plots of the percentage contribution pc(d) to the total entropy as a function of d, for fixed N = 300, n = 
100, 150,200, in each of cases of GS (o = 0) and ES (with o = 30,50). The solid thin curve is for GS (o = 0), 
whereas the bold light and thick curves are respectively for ES with o — 30 and ES with o — 50. 

The variation of the percentage contribution pc(d) with the window width d is shown in Fig. 9 for 
parametric choices of n = 100, 150, 200. In the case of GS, almost the entire entropy is recovered within 
a width of just d = 3 . This again shows that most of the total GS entropy, which obeys the AL, is 
contributed by the interaction region that encompassed the DOF very close to the horizon. In the case 
of ES, for which the AL is violated, it takes d to be as much as 15 — 20, depending on the excitations 
o = 30 — 50, so that the total ES entropy is recovered. Therefore the DOF that are far-away from the 
horizon have a greater contribution to the entropy in the case of ES, than for GS, and such contributions 
increase with increasing excitations o. The location of the horizon (values of n) are also found to affect 
the entropy contributions. In order to estimate of how much the DOF that are far from the horizon affect 
the total entropy in the case of ES, as compared to the case of GS, let us consider the percentage increase 
in entropy for an increment in the interaction region by exactly one lattice point, i.e., 



Apc(d) = pc{d) — pc(d — 1) 



(79) 



This is given by the slope of the pc{d) vs. d plots shown in Fig. 9. In Fig. 10 the variations of Apc(d) 
with (n — d) are shown for n = 100, 150, 200 in the cases of GS and ES (with o = 30, 50). 

For GS, the entropy increases from to about 85% of the total entropy when the first lattice point just 
inside the horizon is included in the interaction region. Inclusion of a second lattice point adds another 
9%, a third lattice point adds 3%, a fourth lattice point adds 1%, and so on. The contributions to the 
entropy by these additional lattice points farther and farther from the horizon decrease rapidly, and by 
the time when the (n/3) th is included the increment in entropy is less than 0.01%. This happens for 
all values of n, i.e., the horizon location does not affect the way by which the DOF contribute to the 
entropy. For ES, however, the inclusion of the first lattice point inside the horizon raises the entropy from 
to about 55 — 75% (depending on n = 100 — 200) for o = 30, and to about 40 — 60% (depending on 
n = 100 — 200) for o = 50. The next successive points add about 9%, 4 — 5%, 3 — 4%, • • • , depending on 
o = 30 — 50 but fairly independent of n, and the corresponding slope is smaller. 
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n = 100 
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n = 150 
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n-d 

Figure 10: Variations of Apc(d) with n — t in the cases of GS (o = 0) and ES (with o = 30, 50), and once again 
for N — 300, n — 100, 150, 200. The solid thin curve is for GS (o = 0), whereas the bold light and thick curves are 
respectively for ES with o = 30 and ES with o = 50. 

From the above results we thus observe that although most of the entropy is contributed by the DOF 
close to the horizon, the DOF that are farther away must also be taken into account for the AL to emerge 
for GS, and the AL plus corrections for ES. With increasing excitations o, the contributions from the 
far-away DOF become more and more significant, as are the corrections to the AL. Thus the AL may 
be looked upon as a consequence of entanglement of the DOF near the horizon, whereas the corrections 
to the AL may be attributed to the contributions to the entropy by the DOF that are far from the 
horizon [15, 16]. 

8 Entanglement entropy of massive scalar field 

In all the analysis presented in this review, we have considered an ideal situation where the scalar field are 
non-interacting and massless. One natural question that arises is what happens if we include interactions 
or if the field is massive? It is, in general, not possible to obtain S cnt non-perturbatively for interacting 
fields and is beyond the scope of this review. In this section, we consider a massive field and show that 
all the results of the previous sections continue to hold. 

The action for a massive scalar field propagating in the background space-time with metric g^ v is 
given by 

S = -±Jd 4 x^ [gT d v y + mV] (80) 

Now, proceeding as before in the case of massless scalar field (see sec. 4), one can obtain the discretized 
Hamiltonian for the case of massive scalar field and show that this Hamiltonian resembles that of a 
iV— coupled HO, Eq. (19), with an interaction matrix i^( m ) whose elements K^ 1 are related to the 
interaction matrix elements Kij [Eq.(20)] for a massless scalar field as 

K \f = Kit + m 2 a 2 . (81) 
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Figure 11: Best fit plots of the relative variation of the total entropy S m for a massive scalar field (in units of the 
total entropy Sq corresponding to a massless scalar field) with the mass m times the lattice spacing a, for fixed 
n = 100, o = 30, in each of cases of GS, ES and MS (Eq/Hi). The corresponding data are shown by asterisks. The 
fits show an exponential damping of the ratio S m /So with mass. 

With this new interaction matrix one then finds the reduced density matrix following the steps 

discussed earlier and finally computes the entanglement entropy S m for the massive scalar field. 

The variation of S m /So [where, So is the entropy due to the massless scalar] with (m x a) is shown 
in Fig. 11 for the cases of GS, ES and MS(Eq/Hi) for fixed parametric values N = 300, n = 100, o = 50. 
The data fit very well with a Gaussian which shows that S m falls off exponentially with respect to Sq as 
the mass increases: 



5,1 



Sq exp 



-ai (ma + a2) 



(82) 



where cvi, «2 and A are the fitting parameters. Depending on the state (GS, ES or MS), the parameter a% 
varies between 2.77 and 3.24, is between 0.077 and 0.001 and the power A is between 2.246 and 2.078. 

2 2 

Therefore, although S m /So approximately scales as e~ m a , from the small variation in the power A one 
finds that the exponential damping is strongest for GS, and gradually slows down as more and more ES 
oscillators are mixed with GS, the damping is slowest for the ES case. 

Even with a fairly high amount of excitation (o = 50) the fitting parameters a%, ot2 and A change very 
little for the different cases GS, MS(Eq/Hi) and ES. Therefore for a particular value of the mass m, the 
relationship between S m and (R/a) practically remains the same as that between Sq and (R/a) in all the 
cases of GS, MS (Eq/Hi) and ES. The analysis and inferences of the previous sections go through for the 
massive scalar field as well, resulting in correction terms as obtained before. 



9 Conclusions 

In the absence of a workable quantum theory of gravity, the best strategy is to slowly build a coherent 
picture and hope to understand — and, in due course, solve — some of the problems of black-hole 
thermodynamics. Thus, it is important to explore all possible avenues. Quantum entanglement, as a 
source of black-hole entropy stands out for its simplicity and generality. The results discussed in this 
review highlight the nontrivial, and somewhat counterintuitive, facets of quantum entanglement and its 
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role as the source of black-hole entropy. More precisely, assuming the modes evolve adiabatically, we have 
shown that: 

• Entanglement leads to generic power-law corrections to the area law 

• The quantum degrees of freedom that lead to 5 BH and subleading corrections are different. 

• It is possible to identify the quantum degrees of freedom that contribute to the area law and the 
subleading corrections. 

• The interactions do not change the form of the entanglement entropy for different quantum states. 

It is important to note that although the analysis presented here is semiclassical, since, the entan- 
glement is a quantum effect and should be present in any theory of quantum gravity. Hence, the results 
presented here do have implications beyond the semiclassical regime. There are some new insights which 
arise in this approach which are worth exploring further: 

• Is there a connection between the entanglement entropy and the Noether charge? On the face of it 
there is no apparent connection: For diffeomorfism invariant theories, like Einstenian gravity, the 
Noether charge is interpreted as the black- hole entropy [13], implying that the higher order space- 
time derivatives contribute to the subleading power-law corrections. In the case of entanglement, 
the subleading contributions arise due to the excited quantum states of the scalar field which exists 
also for black-holes in Einstein gravity. 

However, it has been shown [21] that the classical conserved charge for a nonminimally coupled 
matter fields propagating on the fixed curved background is identical to the Noether charge defined 
in Ref. [13]. This raises the following question: Can the excited states of the quantum scalar field be 
mapped to the non-minimal coupling of the field? If yes, then conserved charge defined in Ref. [21] 
can be related to the Noether charge. This is currently under investigation [22]. 

• Any approach that aspires to explain black hole entropy from fundamental principles must provide 
a natural explaination for the factor 1/4 in the Bekenstein-Hawking entropy. However, in the case of 
entanglement entropy, the proportionality constant in the relation S = 0.3(R/a) 2 for GS obtained in 
ref. [10] differs from the 1/4 in the Bekenstein-Hawking relation [Eq.(2)]. This discrepancy persists 
for MS and ES. A probable reason behind this mismatch is the dependence of the pre- factor on the 
type of the discretization scheme. For example, another discretization scheme, resulting in the NN 
interactions between four or more immediate neighbors, would result in a different pre-factor. Is it 
at all possible to obtain the Bekenstein-Hawking value? 

• In this review, we have assumed that the quantum modes evolve adiabatically thus neglecting the 
contribution of the late-time modes leading to the Hawking particles. What happens if we relax 
this assumption? Will it change the relation (74)? 

The analysis in the previous section suggests that it may not, if we treat the late-time modes 
perturbatively. In Sec. (4), we showed explicitly that the time-dependent Hamiltonian (11) becomes 
a free field Hamiltonian (15) for any Lemaitre time To- At any time t = tq + St, the Hamiltonian 
(11) can be written as 

H BH ~H F + 8H (83) 

where 5H is the small perturbation and includes the interaction. In this case, if (5H)/H F <C 1, the 
relation (74) will not change. However, it will be interesting to investigate the effects at late times 
where the above perturbation expansion fails. 

• Recently, Pleino et al [23] provided analytical proofs of numerical results of Bombelli et al [9] and 
Srednicki [10] and showed that the entropy-area relation do not depend on the shape of the traced 
out volume [23]. It will be interesting to do such an analysis for excited states. 
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• Can a temperature emerge in the entanglement entropy scenario, and if so, then along with the 
current entropy, will it be consistent with the first law of black hole thermodynamics? Are the 
second and third laws of thermodynamics valid for this entropy? Can the entanglement of scalar 
fields help us to understand the evolution or dynamics of black-holes and the information loss 
problem? We hope to report on these in future. 
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A Appendix: Why consider scalar fields? 



In this appendix, we discuss the motivation for considering massless/massive scalar fields for the entangle- 
ment entropy computations, from the perspective of gravitational metric perturbations in asymptotically 
flat spherically symmetric space-times. 

Let us consider the Einstein-Hilbert action with a positive cosmological constant |A|: 



S EH (g) = Ml i j d i x^=g[R-2\h\] . 
Let us decompose the metric g^ v in terms of a background metric g^ v and fluctuations h^ u : 

g^iv = g^v + h[iv ■ 



(84) 



(85) 



Assuming h^ u to be small, and expanding the action keeping only the parts quadratic in h^ u , gives [24] 



S EH {g,h) = -Ml d A x^\g 



I" 



2 7% -T« + \ V„W"fc + ^ h^H" 



where 



(86) 



(87) 



7% = ^(VX + V^-V a V). 

One can easily verify that the above action (86) is invariant under the infinitesimal gauge transformation 
hfj, v — > h^ u + V(^£j,) when the background metric g^ u satisfies the vacuum Einstein's equation in presence 
of the cosmological constant | A| . The gauge arbitrariness can be removed by imposing the harmonic gauge 
condition d~W v = [24]. 

Now by keeping only the first derivatives of h^, the action (86) further reduces to [25]: 



S EH (g,h) = -^Jd 4 x^\g~\ [V a VV a r + |A|V^H- 



(89) 



This corresponds to the action for a massive spin-2 field propagating in the background g^, the mass 
being given in terms of the cosmological constant |A|. 

In the weak field limit can be approximated as a plane-wave perturbation with a particular 
frequency, i.e., 

V = M-h^Mx 11 ) , (90) 

where e^ v is the constant polarization tensor. Consequently, the above action (89) reduces to a form 
which is the same as the action for a massive scalar field ip propagating in the background g^ v : 



S EH (g,h) = -\ Jd 4 x^\g~\ [d a <p&*<p+ \A\ip 2 ] . 



(91) 
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Now, one may further note that in four-dimensional spherically symmetric space-times, the metric 
perturbations are of two kinds — axial and polar [26-28]. The equations of motion of both these pertur- 
bations are scalar in nature and are related to each other by a unitary transformation [26] . The equations 
of motion of the axial perturbations are identical with those of a test, massless scalar field propagating 
in the black-hole background: 

U<p ee ^La M (V=0<rM = . (92) 

Hence, by computing the entanglement entropy of the scalar fields one can obtain the entropy of a 
class of metric perturbations of the background space-time. Of course, such a computation would not 
account for the entropy of all perturbations, because a generic perturbation is a superposition of the plane 
wave modes and the entanglement entropy is a non-linear function of the wave-function. Nevertheless, 
scalar fields are expected to shed important light on the role of entanglement in the AL. 



B Appendix: Hamiltonian of scalar fields in black-hole space-times 



In this appendix, we find the expression for the Hamiltonian of a scalar field propagating in a static 
spherically symmetric space-time and show that for a particular time slicing this Hamiltonian reduces to 
that of a scalar field in flat space-time. 

Let us consider the line-element for a general four-dimensional spherically symmetric space-time: 



ds 2 



df 2 

A(r, dr 2 + + p 2 (r, £) (d9 2 + sin 2 Odtf) , 



(93) 



where A, B, p are continuous, differentiable functions of (r, £). The action for the scalar field p propagating 
in this space-time is given by 



~ J d A x^/^g df.ip d u <p 



-"IE/ 

Im 



drdi 



Vab 



A 



(d T ^ lm y + vabp 2 (d^ lm y + i(i + 1) \i - ^ 



where we have decomposed ip in terms of the real spherical harmonics (Zi m (0, <j))): 

<p{x») = Y l <P, m {T,Z)Z lm (0,<l>). 



(94) 



(95) 



Im 



Following the standard rules, the canonical momenta and Hamiltonian of the field are given by 



TT 8C p2 ft 

d(d T <p lm ) ^AB aTiPl -> 



1 f°° 

Him(r) =2/^ 



n 2 I 



A 



U L + vABp 2 (d£ip tm y + l(l + l)\/-<p 



Im ■ 



(96) 
(97) 



Im 



The canonical variables (ip lm ,U lm ) satisfy the Poisson brackets 

K>,0,n im (r,O} = S(Z - O , W lm (r,0,^ m (r^')} = = {n im (r,0,n im (r,O} • (98) 

In the time-dependent Lemaitre coordinates [25, 29] the metric components of the line-element (93) 
are given by 

A(r,0 = l ; B(t, = T -y^ ; P (r,0=r, (99) 
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where r = r(r, £). The line-element in the Lemaitre coordinates is related to that in the time- independent 
Schwarzschild coordinates, viz., 



ds 2 = -f(r)dt 2 + + r 2 (dfl 2 + sin 2 W^ 2 ) 
/(»") 



/(r = r h ) 



(100) 



by the following transformations [29]: 



r = t± dr 



[1 ~ /(r)]- 1 / 2 
/(r) 



(101) 



Unlike the line-element in Schwarzschild coordinates, the line-element in Lemaitre coordinates is not 
singular at the horizon r^. Moreover, the coordinate £ (or, r) is space(or, time)-like everywhere, whereas 
r(or, t) is space(or, time)-like only for r > r^. 

In the Lemaitre coordinates the general Hamiltonian (97) takes the form 



H lm {r) 



i r 

2Jt 



r 2 \J\- f(r 



-.U 2 + 

Lm 



%j 2 + i(i + i)v^7w^ 



(102) 



which depends explicitly on the Lemaitre time. 

Choosing now a fixed Lemaitre time (r = tq = 0, say), the relations (101) lead to: 

dj _ I 

dr TWW' 



(103) 



If we set d9 = d(j) = 0, then for the fixed Lemaitre time tq it follows that ds 2 = d^ 2 / 'B(tq,£) = dr 2 , i.e., 
the covariant cut-off is \ds\ = dr. Substituting the above relation (103) in the Hamiltonian (102) we get 



H, 



i r 



dr 



" +r 2 (^J 2 + Z(Z + l)yt 



l-/(r 

where the variables (ip lm ,H lm ) satisfy the relation: 

fe>),II ;m (r')} = VT^W)S(r-r'). 
Performing the following canonical transformations 

n. -► rWl-f(r)U, ;ip, -> ^ 
the full Hamiltonian reduces to that of a free scalar field propagating in flat space-time [30] 



lm K 



d ( (fim(r) 



1 2 



+ 



1(1 + 1) 



<Pim( r ) 



(104) 



(105) 



(106) 



(107) 



This happens for any fixed value of the Lemaitre time r, provided the scalar field is traced over either the 
region r G (0, rh] or the region r G [rh, oo). Note that the black-hole singularity can be entirely avoided 
for the latter choice, and for evaluating time-independent quantities such as entropy, it suffices to use the 
above Hamiltonian. 

The approach here differs from that of Ref. [31] where the authors divide the exterior region r > r s 
into two by introducing an hypothetical spherical surface and obtain the entanglement entropy of that 
surface. In contrast, we consider the complete r > r s region and obtain the entropy for the black hole 
horizon. 
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